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ABSTRACT 



. We investigate the delayed, secondary GeV-TeV emission of gamma-ray bursts 

O ■ and its potential to probe the nature of intergalactic magnetic fields. Geometrical 

effects are properly taken into account for the time delay between primary high- 
^ '. energy photons and secondary inverse Compton photons from electron-positron 

pairs, which are produced in 7-7 interactions with background radiation fields 
and deflected by intervening magnetic fields. The time- dependent spectra of the 
delayed emission are evaluated for a wide range of magnetic field strengths and 
redshifts. The typical flux and delay time of secondary photons from bursts at 
z ~ 1 are respectively ~ 10^^ GeV cm~^ s~^ and ~ 10^ s if the field strengths 
are ~ 10^^^ G, as might be the case in intergalactic void regions. We find crucial 
differences between the cases of coherent and tangled magnetic fields, as well as 
dependences on the field coherence length. 

Subject headings: magnetic fields — intergalactic medium — cosmology: the- 
ory — gamma rays: bursts — gamma rays: theory — radiation mechanisms: 
nonthermal 
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Introduction 



Magnetic fields are ubiquitous in the unive rse, being observed no t only in stars but 
also in larger systems, such as spiral galaxies (IBeck fc HoernesI I19961). elliptical galaxies 



( IMoss fc ShukuroyI 119961) .and c lusters of galaxies (JKim et al 



1990 



Carilh fc Tavloi] l2002l ) 



(for a review, see IWidrowl |2002| ) . It has been cla imed recently that even superclusters may 
have associated magnetic fields (IXu et al.l 120061 ) . So far, many observations demonstrate 
that magnetic fields are common in the structured regions. 

On even larger scales, magnetic fields are yet to be observationally determined. A well- 
known technique to probe large-scale magnetic fields utilizes Faraday rotation, i.e. rotation 
of the polarization vector of radiation from background sources by intervening magnetized 
plasma. For example, an upper limit of 6 x 1 0~^^G has be en deriv ed for a uiiiforni magnetic 
field component in the intergalactic medium (IValleelll990l ) (see also lKronberg]|l994l ) . If inter- 
galactic magnetic fields originate from epochs earlier than cosmic recombination, anisotropies 
i n the cosmic microwave background (CMB ) can also provide constraints of nanogauss levels 
(JBarrow et al.l 119971 : lYamazaki et al.ll2006l ). 



Various theoretical possibilities exist for the origin of magnetic fields on very large 
scales. During the inflationary epoch, magnetic fields on cosmological s cales can be gener- 
ated i f conformal i nvariance of electromagnetic in t eract ions is violated (JTurner fc Widrow 
19881 : iRatral Il992l : iBamba &: Sasaki! 120071 : iBambal 120071 ). Large-scale magnetic fields can 
also arise through the evolution of priniordial density fluctuations bef ore cosmic recombina- 



tion (Takahas 



2007al : iHogan 



li et a. 



2000 



2OO5I: llchiki et al.ll2006bl. I2OO7I: iKobavashi et al 



Berezhiani &: Dolgov 



2004 



Matarrese et al 



boosh " 



20071 : iTakahashi et al. 



Although the field 

amplitudes are not expected to be large, < 10~^° G, their generation is inevitable, being in- 
duced by the observed density fluctuations. Even after recombination, it has been proposed 
that magnetic fields of ~ 10~^°-10~^^ G can emerge at cosmic rei onization fronts aro und 
redshift 2; ~ 15, throu gh either the Bierm ann battery mechanism (JGnedin et al.l l2000l ) or 
radiation drag effects (JLanger et al.ll2005l ). Such large-scale magnetic fields originating at 
high z may remain very weak in some intergalactic regions at lower 2;, well below the upper 
limits from Faraday rot ation, particularly in void regions where the activity of astrophysical 



objects is minimal (e.g. iBertone et al.ll2006l ) 



Plagal (I1995I ) proposed a promising method to observationally probe such weak inter- 



galactic magnetic fields, relying on the delayed, secondary high-energy emission from gamma- 
ray bursts (GRBs). If primary photons of the GRB prompt emission have sufficiently high 
energies, they can interact with photons of the cosmic infr ared background (G I B) or CMB to 



creat e electron-positron pairs far away from the GRB (e.g. IStecker et al.lll992l : iKneiske et al. 



20041 ). These charged particles are then deflected by intergalactic magnetic fields before 
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emitting secondary GeV-TeV gamnia-rays via inverse Compton (IC) upscattering of the 
CMB, which reach the observer with a time delay relative to the primary gamma-rays. This 
delayed emission depends on the propertie s of the intervening ra agnetic fields and hence 
constitute a valuable probe of its nature. iRazzaque et al.l (120041 . hereafter RMZ04) have 
discussed the time-dependent spectra of suc h delayed emiss i on components based on sim 



pie est i mates of the delay t imescale (see alsolDai fc Lull2002l : IWang et al.ll2004j : llchiki et al. 



2006al : iMurase et al.l 120071 : ICasanova et al.l 120071 . for similar and related discussions) . Al- 
thou gh observatio r ial information on GeV-TeV emission from GRBs has been limited so far 



[e.g. lDingusll200ll : lAlbert et al.l 120071 ). significant advances are expected soon from the new 
generation of satellite- and ground-based facilities that have recently begun operation or will 
come online in the near future. 

Here we reconsider this problem, with the aim of achieving a more physically consistent 
picture of the delayed emission, as well as clarifying the dependence on the magnetic field 
properties in greater detail. First, we develop a new formalism that properly accounts 
for geometrical effects in the delayed radiation, with important consequences for the time- 
dependent spectra. Second, whereas previous studies treated only coherent magnetic fields, 
we discuss cases of both coherent and tangled magnetic fields, and the resulting crucial 
differences. We also explore the dependenc e on redshift through the evolution of the GIB, 
adopting the GIB models of iKneiske et al.l (|2002| . 120041 ). The observational prospects are 
briefly addressed. 

In §2, we describe our new formulation. The resultant spectra of the delayed emission 
are presented in §3. We provide a discussion in §4, and conclude with a summary in §5. 



Formulation 



2.1. Preliminaries 



Our new formalism for the delayed gamma-r ay eniission from GRBs ern ploys a simple 
description of the GRB prompt emission (see, e.g. |Piranll2005l : lMeszarosll2006l . for reviews on 
GRBs), together with an adaptati on of a treatment for dela yed radiation due to scattering 
effects in the ambient medium (e.g. lSazonov fc Sunyaevll2003l ). Following RMZ04, we assume 
a simple power-law spectrum of photon index a for the high-energy part of the prompt 



emission, 



dE^dt 



(a 



1)L 



7, ISO 



(t) 



^^DlEl^^ 



E^ 



E. 



'7>pk 



{E 



7,pk 



< E^ < E, 



cut/ 
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where L^^iso{t) is the isotropic-equivalent gamma-ray luminosity in the observer frame, -E'7,pk 
is the observed spectral peak energy, typically a few hundred keV, and E^nt is the maximum 
cutoff energy. We allow for time-dependence in the luminosity, although the spectral shape 
is assumed to be constant for simplicity. The luminosity distance is given by 

Dl{z) = ^^^^ rdz'[nr{i + zT + ^m{^ + zT + nAr'^' (2) 

-no Jo 

where Qr, ^m and ^a are the densities of relativistic particles, non-relativistic matter and 
cosmological constant in units of the critical density, respectively, and Hq is the Hubble 
constant. 

The primary high energy photons >100 GeV interact with GIB and CMB photons to 
create high energy electron-positron pairs (e^). These e^ are responsible for the delayed 
emission of gamma-rays through IC upscattering of CMB photons. By assuming that each 
electron or positron (hereafter simply "electron" ) share half the energy of the incident pri- 
mary photon, me7e = E.y/2, the flux of e^ corresponding to time t = tcRB during the burst 
duration can be written as 



dN^ = "-^-''r'" "!-. -£ Il - exp(-r(27.m,))| . (3) 

d-fedtcRB ^T^Dl (2me)" ^^^,pk 

where t{E^) is the optical depth to pair production with the CIB and CMB for gamma-rays 
of energy E^. 



While the number density of CMB photons is known preci sely, that of th e CIB 



arising 



from stellar and dust-reprocessed emission is still uncertain. iKneiske et al.l (120021 . |200J) 
have developed semi-empirical, backward evolution models of the CIB based on observed 
data from galaxy surveys. We adopt their "high-stellar- UV" model which gives the best 
description of the observed proximity effect in high- 2; quasar absorption line systems, even 
though their "best-fit" model also leads to essentially the same end results. In Fig. [H we 
show the optical depth to pair production employed in this paper. Note that the two peaks 
in the CIB curves each correspond to the contributions from stars and dust. 



Once one obtains the e^ energy distribution, the spectrum of the delayed emission can 
be calculated as 

j2 a 

^ delayed 



d^Ndel^yed f , dN^d^N^c 



d'j^—-^ ^^ . (4) 

dtdE^ J ^ d'^f. dtdE^ 

Here dNe/d'-je is the total time-integrated flux of e^ responsible for the delayed emission 
observed at time tobs after the burst, and 

d^Nic 3aT /• , , , fix) 

'dME, - A^'J ^^7,CMBncMB(e,,cMB)^-^ , (5) 
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Fig. 1. — Optical depths due to the cosmic infrared background (IR; red-dashed) and cosmic 
microwave background (CMB; black-solid) at redshifts z = 0.1 to 4.0. 



represents the IC power from a single electron in the Thomson regime, where ecMB and 
n(ecMB) are the energy and nur nber density per unit energy of CMB photons, respectively. 
The function /(x) is defined as (JBlumenthal &: Gouldlll970l ) 



/(^) 



2a; In X + X -|- 1 — 2x^ 



X 



47eeCMB 



(6) 

(7) 



At this point we need to relate dNe^o/d-fe in Eq.© with dNe/d-fe in Eq.(jl]). RMZ04 
use a simple approximation and assume dNe/d'je = {tic/ ^tB)dNefi/d'je, implying that the 
observer receives continuously for a typical delay time Ate all photons that were emitted by 
e^ during the IC cooling time tic- However, a more consistent treatment should consider 
the proper geometrical configuration of the delayed emission. The primary photons travel 
distances of several Mpc or more depending on their energy before creating e^, which then 
emit secondary photons while propagating outwards, subject to deflection by intervening 
magnetic fields. In the following we formulate a relation between dN^fi/d'ye and dN^/d'ye, 
taking into account these crucial geometrical effects. 

Note that the delayed secondary emission of Eq.(j4]) itself suffers 7-7 absorption with 
the CIB and CMB according to t{E^) before reaching the observer, which we take into 
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account. However, we neglect additional cascading effects that involve further generations 
of e^ production and IC emission. The validity of this approximation is addressed later in 



2.2. Geometry of delayed emission 

Let us suppose a GRB at redshift z radiates the spectrum of Eq.([T]). The geometrical 
configuration is depicted in Fig. [2l A high-energy primary photon emitted in the direction 
of Q with respect to the line of sight travels a distance of Air before interacting with a GIB 
photon to create an e^. The 7-7 mean free path in the GIB is 



^» = osw;; == 2° '''p^ ( 



0.26(Jt^ir V0.1cm~3 

where gt is the Thomson cross section, and rii^ is the number density of GIB photons 
appropriate for given E^ and 2;, the numerical value corresponding to E^ ^ 10 TeV and 
z = 1. The electron then upscatters ambient GMB photons with IG mean free path iic 
multiple times to produce secondary emission until it loses most of its initial energy after 
propagating an IG cooling length Ajc, which can be written 



1 



!^IC 



^lkpc(l + z)"^ (9) 

O'T'ncMB 

Aic = , ^""^ - 35 kpc (tttItt) V + ^)-^ (10) 

47eCrrpcMB \W TeVy 

where ncMB and pcmb are the number and energy density of GMB photons, respectively. 
During propagation, the electrons are deflected gradually by successive IG scatterings as 
well as by the ambient magnetic fields. The deflections can be modeled as a random walk in 
angle, and the probability that an electron is deflected by an angle 6' from its initial direction 
0is 

m 0') = -^ ^^- exp (- ^^'~fjA , (11) 

which is normalized so that j PdQ' = 1. Here a{6) and (6) are respectively the variance 
and expectation value of the deflection angle in the 6 direction, which encode important 
information about the magnetic flelds such as their amplitude and coherence length. They 
are discussed below separately for the cases of tangled and coherent magnetic flelds. 

The secondary photons will be observed if their parent electrons have been deflected by 
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constant time delay surface ' 







A^^-^"' ,' 



low energy photon 



Fig. 2. — Schematic picture of the geometrical configuration. Left: Overall geometry. The 
angles 6'min and 6'max delineate the range of integration in Eq. flT2|) . Right: Expanded view 
of the left panel. "Effective constant delay surface" is the surface from which the delayed 
photons reach the observer at the same time, including the delay during propagation through 
intervening magnetic fields. See main text for details. 



an angle ^, and we can express dN^/d'^e as 



dN, 

die 



dt 



GRB 



(*obsi*GRB) 



'max (,robs ) 



6min(*obs>tGRB) 
Smax(tobs>tGRB) 



de / dn'5{e - e')p{e\ 



djedt, 



GRB 



dt 



GRB 



de 



exp 



{9-{9)f\ dKfl 



2a\e) / d-iJtGKB 



(12) 



^mi„(tobs,tGRB) y2'Ka{e) 

where the directional requirement is represented by the delta function 5{9 — 9'), and the 
range of integration [6'mm(^obs,^GRB), ^max(i^obs,^GRB)] depends on the observation time (see 
Fig. [2]). In the second equality of Eq. 0121) . the integration over ^ dVt' = ^ dcj)' s\.n9'd9' 
was performed. The observed delayed emission can then be expressed as an integration over 

dNe/d-ie, 



d^N, 



delayed 



dtohsdE^ 



d-f, 



dNf. 3o"t d (r) 



'^e7,CMB^CMB(e7,CMB) 



e7,CMB 



(13) 



Since all photons emitted frona the constant delay surface reach the observer at the same 
time (ISazonov fc Sunyaevll2003l ). here we have replaced c in Eq.([5]) with d (r) /dtohs, where 
(r) is the rectilinear distance between the locations of the electron at production and at tobs 

(Fig. m- 



2.2.1. Tangled magnetic fields 



If the ambient magnetic fields are tangled, that is, their coherence length is much smaller 
than the electron cooling length, electrons undergo a random walk in angle not only by IC 



scattering but also by magnetic deflection. For such fields, 

Ao) = {om) + (efM) + ^2 ^ (14) 

(^) = 0. ' (15) 

The last term in Eg. (1141) describes the effect of angular spreading by the initial pair produc- 
tion interaction. The variance of the deflection angle due to magnetic deflection {0%{0)) and 
that due to angular spreading by IC scattering {61q{6)) are given by 





1 


i{oim - 


Tc R{0) 


'{Ohm - 


ecMB / -Kl^yj 
mec2 V 3£ic 



(16) 

(17) 

where r^ = jgrnec/eB is the Larmor radius of the electron in a magnetic fleld with amplitude 
B, re is the coherence length of the magnetic fleld, and R{0) is the distance the electron would 
travel in the absence of any deflection (Fig. [2]). The relation between (r) and R{0) is 

R-{r) = ^R{<t>l) + ^R{4>!c) , (18) 

where (^b and tb are respectively the deflection angle and "optical dept h" due to magnetic 



deflection, and 0ic and ric are those due to IC scattering (see Eq.(9) in lAlcock fc Hatchett 



19781 ). With the approximations 



0B ~ — , Tb^ — , (19) 

0IC ~ - , ric ~ ^ , (20) 
Eq. (IT8l) becomes, 

(r) = i^-^e^ (21) 

0' = (^l(^)> + (^i'c(^)>- (22) 
Meanwhile, tobs is related to R and (r) as 
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Using these relations and the fact that 0^1, R{6) can be expressed in terms of tobs and 6 
as 



R "^ 



92 COS ^ 



02 fe^\\e^cosefUs-tGRB o\ 

T + vUJ +^:^v i+z "2 ^^ 



(24) 



2 ftohs — ^GRB d . 



e"^ \ i + z 2 

For the range of integration in Eg. (1121) . the condition (r) = /? = in Eq. (l23|) together 
with 6'max ^ 1 gives 

V Air(1 + z) 



On the other hand, taking R = Aic in Eq. fl25l) leads to 



/2(tobs - ^gRb) /^„x 

where Atot = Aic + Air. The expression for d (r) /dtohs can be obtained as follows. Differen- 
tiating Eqs. fl23|) and fl2Tl) with respect to tobs, 

1 di? d (r) 



1 + ^ dtobs C?t 



■cos^, (28) 



obs "'''obs 



''^ - ^-ie^V'£««fi + ie^^l«. (29) 



dtobs V 6 / (itobs V 6 / (itobs 

From these equations one obtains 

d{r) _ 2c 

dtZ~ {l + z) [6^ + 673] • 

In summary, with the above formulation, the factor 



(30) 



, ''"' 2 (a-(fl))^ d (r) 

S= / rf^_J e"^;^^(^^ (31) 



should be integrated, instead of c(tic/AtB) as in RMZ04. 
In what follows, for simplicity we take in Eq. ([3j) 

-^7,iso('^GRB) = -^7,iso^(^GRB)'A^GRB , (32) 

implying that all primary photons are emitted at ^grb = 0, and L^^iso is now the luminosity 
averaged over the duration of the prompt emission AtcRB? a good approximation when 

^obs ^ AtcRB- 



-10- 



2.2.2. Coherent magnetic fields 

If the ambient magnetic fields are coherent, or their coherence length is much larger than 
the electron cooling length, they do not contribute to random walks in angle. We assume 
then that the magnetic deflection can be described by a single, small angle scattering, so 
that Eqs. (ll4lll5p are replaced with the approximations 



a\9) 



ic 



W> + 



I'e 



(^) = - 



R 



as well as Eq. (l5Ul) with 



d{r) 
dt, 



2c 



obs 



[l + z)[9^ + {9U0))/S] ' 



(33) 
(34) 

(35) 



[e.g. iSazonov fc Sunyaevll2003l ). The resulting differences in the delayed emission from the 
tangled field case will be crucial for probing the nature of the magnetic fields, as discussed 
below. 



3. Results - spectra of delayed emission 

The time- dependent spectra of the delayed secondary emission can be calculated with 
the formalism developed in the previous section. Because our main interest is in probing the 
state of the intervening magnetic fields, first we fix the GRB model parameters to a fiducial 

set. 



(L^,i,o,a,^^,pk,^cut,AtGRB) = (10^'erg, 2.2, 500 keV, 10 TeV, 50 s) 



(36) 



most of which are typic ally observed values for the prompt emission ( JFishman &: Meegan 



19951 : iPreece et al.ll2000l ). The only exception is -Ecut? for which the current observational 



constraints are quite poor in the GeV-TeV bands due to the limited sensitivity and/or field 
of view of previous generation instruments, as well as the strong effects of intergalactic 7-7 
absorption for typical GRB redshifts (Fig. [T]). However, in the widely discussed, internal 
shock model of GRB prompt emission, one c an expect that at least some bursts have spectra 



extending to TeV energ i es and above (e.g. IZhang fc Meszarod |2004J : iRazzaque et al.l 12004 



Dermer fc AtoyanI l2006l : lAsano &: Inoud 120071 ) . which we assume to be the case here. The 
consequences of varying the GRB parameters are discussed in ^ 

For the cases of random magnetic fields, we set Tc = 100 pc as a reference value for the 
coherence length. Such values may be related t o the smallest scale of fields generated during 
cosmic reionization by radiation drag effects (JLanger et al.ll2005l ). or to those observed in 
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clusters of galaxies (~ kpc; IVogt fc EnfilinI ( 20051 )). A further, practical reason is that Tc 
should be much smaller than the IC cooling length (Eq JTOj) in order for our random walk 
formulation to be valid. 
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Fig. 3. — Spectra of delayed emission for a fiducial GRB at 2; = 1 and tobs =100, 10^, and 
10^ s, for tangled magnetic fields with B = 10^^^ G and Vc = 100 pc. Thin (blue) dashed 
and solid lines are the spectra of the prompt emission before and after attenuation by the 
GIB, respectively. 

Fig. [3] shows the time- dependent spectra of the delayed emission for a fiducial GRB at 
2; = 1 in the case of tangled magnetic fields with B = 10~^^ G and re = 100 pc. The spectra 
generally consist of a hard power-law part at low energies with photon indices < 2, typical 
of IG emission from e^ injected at high energies, and a steeply falling part at high energies 
(> 100 GeV for z = 1) due to 7-7 attenuation by the GIB. The flux at high energy decays 
with time much more rapidly than at low energy, resulting from the shorter IG cooling time 
as well as the shorter delay time due to magnetic deflections for the corresponding higher 
energy e^. These properties are distinctive characteristics of the delayed secondary emission. 

The delayed emission spectra for different magnetic field strengths are displayed in Fig. 
m for the case of tangled fields with Vc = 100 pc. Although changes in the intervening 
magnetic fields can significantly affect the time evolution of the delayed emission, the total 
fiuence will remain the same. Stronger fields induce larger delays so that the emission has 
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Fig. 4. — Spectra of delayed emission for a fiducial GRB at 2; = 1 and tobs 
10^ s, for tangled magnetic fields with r^ = 100 pc and B — t^^^'^ 



=100, 10^ and 



10-^■^ 10-l^ and 10-19 G. 



lower flux at early times but lasts longer, whereas weaker flelds result in higher initial flux 
that decays quickly (see also RMZ04). If the flelds are sufficiently weak, B < lO^^'^G, the 
delay due to angular spreading by pair creation and IC scattering dominates over magnetic 
deflections, and the behavior becomes independent of B. Note that B and r^ for tangled 
flelds enter in our formulation only through the combination B^frl (Eq. [T6l) . so varying Tc 
leads to results analogous to varying B. 

Fig. [5] is a comparison between the cases of tangled and coherent magnetic flelds with the 
same fleld strength. We see that the spectral evolution is markedly distinct between the two 
cases, owing to different dependences on 7e expected for the mean deflection angle: 



7e for tangled flelds, as opposed to Qb oc 7"^ for coherent flelds. This offers interesting 
prospects for observationally probing not only the fleld strength but also the coherence 
length, providing valuable insight into the physical nature and origin of the intergalactic 
magnetic flelds. 

In Fig. El we depict how the spectra at flxed tobs varies with redshift. Besides the 
obvious dependence of the flux on luminosity distance, the spectral shape of the steep, high- 
energy part changes considerably due to large differences in the amount of 7-7 attenuation. 
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Fig. 5. — Spectra of delayed emission for a fiducial GRB at z = 1 and different tobs as 
indicated, for coherent magnetic fields with B = 10~^^ G (red dashed lines), and tangled 
magnetic fields with B = 10~^^ G and Vc = 100 pc (black solid lines). 

Also shown in Fig. [6] are th e sensitivies for GLAS T [^, and for the Cherenkov telescopes 
MAGIC, H.E.S.S. and CTA JPe Angelis et allboOTl ). Note that these are actually limiting 
fluxes for steady sources with the given integration times, but should serve as guides to the 
detectability, considering that the delayed emission decreases monotonically with time and 
the fluxes at earlier times are always larger. The observational prospects are interesting for 
GLAST and current Cherenkov telescopes for bursts at z < 1, and even a.t z > 1 for the 
future CTA project. 

Finally, let us comment on the difference of our results with those of RMZ04. They 
evaluate the typical delay timescale due to magnetic deflection as Ate = ^ic6'b/2, where 9b = 
Aic/tl is the deflection angle in a coherent magnetic field. More appropriately, however, the 
typical delay timescale should be given by Ate = {tic+^m/cj^B/'^y because the primary high- 
energy photons propagate a significant distance of Air before creating e^ (Fig. [2]). Neglecting 
this crucial geometrical aspect will underestimate the effects of magnetic deflection on the 



^http://glast. gsfc.nasa.gov 
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Fig. 6. — Spectra of delayed emission for a fiducial GRB at tobs 



10^ sand z = 0.1, 0.5, 



and 1.0, for tangled magnetic fields with B = 10 ^^ G and Tc = 100 pc. Overlayed are 5a 
sensitivities for GLAST, MAGIC, H.E.S.S. and CTA with integration times of 10^ s. 

delayed secondary emission. 

We have calculated the delayed emission spectra in a modified version of the RMZ04 
formulation by replacing Ate = ticO% with Ate = (^ic + Mr/c)9%/2. The redshift evolution 
of the CMB and GIB are also considered, which we believe were not properly accounted for in 
RMZ04. The results are compared with those of our full calculation for coherent magnetic 
fields in Fig. [71 showing that they are broadly consistent with each other, despite some 
remaining differences at low energies and at early times. 



Discussion 



As demonstrated above, the spectral evolution of the delayed secondary emission from 
GRBs depends not only on the strength of the intervening magnetic fields, but also strongly 
on whether the fields are tangled or coherent, and for the former case, on the coherence 
length as well. In general, magnetic fields comprise components with a range of coherence 
lengths, characterized by its power spectrum. When the spectrum can be approximated by 
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Fig. 7. — Spectra of delayed emission according to a modified version of RMZ04 (black solid 
lines) compared with our full calculation (red dashed lines), for a fiducial GRB at z = 1 and 
tobs = 10^ 10^, 10^ for coherent magnetic fields with B = IQ-^^ G. Note that the RMZ04 
curves for tobs = 10^ and 10"^ are indistinguishable. 



a simple power-law, the type of field component that determines the magnetic deflection 
angle depends on the spectral index ub- Denoting the spectrum as function of scale by 
B(X) = y/{B^)^ oc A"^, the deflection angle by coherent components can be expressed 

as 6'coh oc i?(Aic)Aic, while that by tangled components is -w (^rand) oc B{rc)Vr\Jvic (for 

re < Aic)- From these relations, 6'coh/^rand = (Aic/''"c)"^^ , so that tangled flelds dominate 
when ub < —1/2. The time-dependent spectra discussed in ^ for tangled and coherent 
flelds should correspond to cases with ub ^ —1/2 and ub ^ —1/2, respectively. 

Let us now briefly discuss the observational prospects. If magnetic flelds in intergalactic 
regions of 7-7 e^ production are too strong, the time delay can become so large that the flux 
of the delayed emission becomes too low to be detectable by GeV-TeV instruments. This 
implies an upper bound to the magnetic flelds strengths that can be probed, depending on 
the brightness of the GRB and the sensitivity of available observational facilities. On the 
other hand, as mentioned in §31 a lower bound of i? ~ 10"^'' G is set by the requirement that 
the time delay is determined by magnetic deflection rather than by angular spreading from 
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IC scattering. Since the higher energy delayed photons generally come from higher energy e^ 
that experience smaller deflections, they may be more suitable for probing stronger magnetic 
fields in comparison to lower energy delayed photons. Conversely, the lower energy delayed 
photons may have better chances to probe weaker fields. Such issues related to o bservability 
will be addressed in more detail in a future publication (JTakahashi et al.ll2007bl ). 



We emphasize that the geometry of the processes involved in delayed secondary emis- 
sion is naturally configured to probe only the magnetic fields in intergalactic locations such 
as voids where they are likely to be weak, without being affected by stronger fields existing 
around the structured regions of the universe. First, the primary photons travel distances 
of Air ~ 20 Mpc (niR/0.1cm^)~^, the 7-7 mean free path in the CIB, depending moderately 
on their energy for E^ ~ 1-10 TeV. This should be far outside the GRB's host galaxy, and 
presumably out into intergalactic void regions away from large-scale structure. Furthemore, 
once they turn into e^ and become subject to magnetic deflections, they cool relatively 
quickly and stop radiating on scales of Aic ~ 35 kpc (-Ee/lOTeV)"^, so that the associated 
secondary emission is sensitive only to the fields within the void regions. Finally, the sec- 
ondary photons will be indifferent to any foreground extragalactic or Galactic magnetic fields 
as long as they do not induce further e^ generation (see below). These features are of great 
advantage over Faraday rotation methods, which can only provide measures integrated along 
the line of sight and are always contaminated by the Galactic contribution. 

It remains to be seen, however, whether the field strengths in such intergalactic voids 
are in the observationally favorable range of i? ~ 10~^°-10~^^ G as discussed here. It could 



plausibly be the case depending on their physical origin at high re dshift (e.g. iGnedin et al. 



2OOOI : iLanger et al.l l2005l : llchiki et al.l 120071 : iTakahashi et al.l l2007al ) , as long as they are not 
significantly polluted later by magnetiza tion from astrophysical so urces such as galactic winds 
(JBertone et al.ll2006l ) or radio galaxies (JFurlanetto fc Loebll200ll ). We stress that such weak 
intergalactic fields would be very difficult to constrain observationally other than through 
the delayed GeV-TeV emission. 

In relation to the observability, we caution that the GRB itself can possess long-lasting 
after g low emission in the GeV - TeV bands (e.g. iMeszaros fc Reed 1 19941 : iBottcher &: Dermei 



19981 : IZhang &: Meszarosll200ll : llnoue et al.ll2003l ). which may potentially mask the delayed 
secondary emission of our interest here. However, the delayed emission should be distinguish- 
able from the afterglow through clear differences in their time evolution. Fig. [8] displays the 
light curves of the delayed emission for the case of Fig. [3], which decay exponentially at late 
times corresponding to the typical delay time Ate = (^ic + '^ir/c)6'^/2, with a characteristic 
energy dependence. In contrast, power-law declines are generally expected for afterglow light 
curves until late times, apart from possible breaks due to jet expansion, etc. If the GeV-TeV 
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afterglow is domi nated by IC emission fro m the forward shock, which could be the most 
common case (e.g. lZhang &: Meszarosll200ll ). simultaneous measurements of the synchrotron 



emission in the radio to X-ray bands will also provide tight constraints on the behavior of 
the high energy afterglow component. 
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Fig. 8. — Light curves of the delayed emission at E^ = 1, 5, 10, and 50 GeV, for a fiducial 
GRB at 2; = 1, for tangled magnetic fields with B = 10^^^ G and Vc = 100 pc. 

A further caveat regards uncertainties in the spectrum of the primary emission. Fig. i^ 
shows how the delayed emission spectra depends on a and E^nt of the GRB prompt emission 
as described by Eq.([T]). Varying a and/or i^cut ^ 10 TeV have large effects, as it changes 
the amount of primary photons above the threshold for 7-7 interactions with the GIB. It 
appears to be insensitive to changes of Ecut ^10 TeV, but this is actually an artifact of 
our approximation where we neglected the effect of cascading beyond the second generation 
of e^ production, which could in reality be important for high -Ecut- Such multi-generation 
cascades entails a complicated geometry that cannot be treated analytically and requires 
numerical studies. We remark, however, that i^cut is probably not much greater than ~ 10 
TeV, or else it would imply uncomfortably lar ge bulk Lorentz factors for the GRB outflow, 
independent of the emission mechanism (e.g. iLithwick fc Saril 1200 ll : iBaringI 120061 ). Note 
also that internal shock models generally predict separate, high-en ergy spectral components 



above a simple powe r-law extrapolation from low energies (e.g. IZhang fc Meszarod 12004 



Asano &: Inoud 120071 ) which would enhance the delayed emission, so the characterization 
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of Eq.([T]) may be too simplistic. More detailed and quan titative studies of the problem 
including such realistic aspects will be presented elsewhere (JTakahashi et al.ll2007bl ). 
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Fig. 9. — Spectra of delayed emission for a GRB at z = 1 and tobs = 10^ s in tangled 
magnetic fields with B = 10"^*^ G and Vc = 100 pc for different a and Ecut as indicated, 
other GRB parameters remaining fiducial (Eq.( !36l) ). 

Finally, we comment that the present approach to probing intergalacti c magnetic fields 



may also be ap plicable in some instances to flaring emission from blazars (jDai et al.l 12002 



Fan et al.l 120041 ). Recent observations of a nearby blazar by the H.E.S.S. collaboration have 
revealed e xtraordinary fiaring ac tivity on minute timescales with spectra extending to sev- 
eral TeV (jAharonian et al.l 120071 ). Searching for delayed GeV- TeV emission from such ob- 
jects may offer an important probe of weak intergalactic magnetic fields in foreground voids. 
Although their nontrivial light curves pose complications, blazars have some crucial advan- 
tages over GRBs. Besides the much larger photon statistics achievable, the intrinsic primary 
emission before 7-7 attenuation can be much better cons trained by virtue of simul taneous 
observations of the lower energy, synchrotron component (ICoppi fc AharonianI 119991 1. 



Observational tests of intergalactic magnetic fields through GeV- TeV gamma-ray as- 
tronomy should be forthcoming with the anticipated launch of GLAST, and ongoing devel- 
opments with ground-based instruments such as MAGIC, H.E.S.S., VERITAS, CANGAROO 
III, MILAGRO, etc. The typical fluxes discussed above may be particularly interesting for 
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the upgraded MAGIC II and H.E.S.S. II telescopes with their lower threshold energies, not 
to mention future projects such as CTA and 5@5 (Fig. [6]). In-depth discussions of the 
detectability with such facilities is deferred to future work. 



5. Summary 

With a new formulation that accounts for the proper geometry of the problem, we 
have explored in some detail the delayed, secondary GeV-TeV emission from GRBs and its 
potential to probe key properties of weak intergalactic magnetic fields. The delayed emission 
was found to possess unique characteristics that depend strongly on the amphtude as well 
as coherence length of the magnetic fields. In the case of intergalactic fields with strength 
B ~ 10~^® G and coherence length Tc ~ 100 pc, the fiux expected from a typical GRB 
at 2 ~ 1 is ~ 10~^ GeV cm~^ s~^ at ~ 10^ s after the burst, in the range of interest for 
current and near future observational instruments. GeV-TeV observations of GRBs may 
soon constrain the weakest intergalactic magnetic fields, which could be difficult to do in 
any other way, offering valuable insight into the physics of the intergalactic medium at high 
redshift, and possibly into cosmological processes in the early universe. 

We are very grateful to Tanja Kneiske for making her GIB models available on her 
website. KT thanks K. Murase and S. Nagataki for helpful suggestions and useful discussions. 
SI thanks P. Goppi for valuable comments. KT and KI are supported by a Grant-in-Aid for 
the Japan Society for the Promotion of Science Fellowship. 
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